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(Dated: January 16, 2008) 

Abstract 

Using the CLEO-c detector at the Cornell Electron Storage Ring, we have measured inclusive 
and exclusive cross sections for the production of D+, and Df mesons in e~^e~ annihilations at 
thirteen center-of-mass energies between 3.97 and 4.26 GeV. Exclusive cross sections are presented 
for final states consisting of two charm mesons {DD, D*D, D*D*, DfDj, Dl~^Dj, and D*'^D*~) 
and for processes in which the charm-meson pair is accompanied by a pion. No enhancement in 
any final state is observed at the energy of the F(4260). 
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I. INTRODUCTION 



Hadron production in electron-positron annihilations just above cc threshold has been a 
subject of mystery and little intensive study for more than three decades since the discovery 
of charm. Recent developments, like the observation of the y(4260) reported by the BaBar 
collaboration [1] and subsequently confirmed by CLEO-c [2] and Belle [3], underscore our 
incomplete understanding and demonstrate the potential for discovery of new states, such as 
hybrids and glueballs. It is also clear that precise measurements of charm-meson properties 
will shed light on higher-energy investigations of 6-flavored particles and new states that 
might decay into b. Charm decays also offer unique opportunities to test the validity and 
guide the development of theoretical tools, like lattice QCD, that are needed to interpret 
measurements of the CKM quark-mixing parameters [1]. Any comprehensive program of 
precise charm-decay measurements demands a detailed understanding of charm production. 

Past studies of hadron production in the charm-threshold region have been largely mea- 
surements of the cross-section ratio R{s) = a{e~^e~ —>■ hadrons)/(T(e"'"e~ —>■ fi'^fJ'~) over this 
energy range that have been made by many experiments Recent measurements with the 
Beijing Spectrometer (BES) [6] near charm threshold are especially noteworthy. There is a 
rich structure in this energy region, reflecting the production of cc resonances and the cross- 
ing of thresholds for specific charm-meson final states. Interesting features in the hadronic 
cross section between 3.9 and 4.2 GeV include a large enhancement at the threshold for 
D*D* production (4.02 GeV) and a fairly large plateau that begins at D*~^D~ threshold 
(4.08 GeV). While there is considerable theoretical interest [3 El El ED] , there has been little 
experimental information about the composition of these enhancements. 

In this paper we describe measurements of charm-meson production in e"'"e~ annihilations 
at thirteen center-of-mass energies between 3970 and 4260 MeV. These studies were carried 
out with the CLEO-c detector at the Cornell Electron Storage Ring (CESR) [II] in 2005-6. 
(Throughtout this paper use of any particular mode implies use of the charge-conjugate 
mode as well.) The principal objective of the CLEO-c energy scan was to determine the 
optimal running point for studies of D^-meson decays. The same data sample has been 
used to confirm the direct production of y(4260) in e^e~ annihilations and to demonstrate 
y(4260) decays to final states in addition to n^n^ J/ip [2\. Specific results presented in this 
paper include cross-section measurements for exclusive final states with D~^, and Df 
mesons and inclusive measurements of the total charm-production cross section and R{s). 

II. DATA SAMPLE AND DETECTOR 

The data sample for this analysis was collected with the CLEO-c detector. Both the 
fast-feedback analysis carried out as data were collected and the detailed analysis reported 
here are extensions of techniques developed for charm-meson studies at the V'(3770) [12]. 

An initial energy scan, conducted during August-October, 2005, consisted of twelve en- 
ergy points between 3970 and 4260 MeV, with a total integrated luminosity of 60.0 pb~^. 
The scan was designed to provide cross-section measurements at each energy for all acces- 
sible final states consisting of a pair of charm mesons. At the highest energy point these 
include DD, D*D, D*D*, D+D;, D*+D;, and D*+D*-, where the first three include both 
charged and neutral mesons. A follow-up run beginning early in 2006 provided a larger 
sample of 178.9 pb"^ at 4170 MeV, not one of the original scan points, that proved essential 
in understanding the composition of charm production throughout this energy region. 
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The center-of-mass energies and integrated luminosities for the thirteen subsamples are 
hsted in Table [ij Integrated luminosity is determined by measuring the processes e~^e~ — > 
e"'"e~, f^'^fJ,~, and 77 ^13j, which are used because their cross sections are precisely determined 
by QED. Each of the three final states relies on different components of the detector, with 
different systematic effects. The three individual results are combined using a weighted 
average to obtain the luminosity used for this analysis. 

CLEO-c is a general-purpose magnetic spectrometer with most components inherited 
from the CLEO III detector [Mj, which was constructed primarily to study B decays at the 
T{4S). Its cylindrical charged-particle tracking system covers 93% of the full An solid angle 
and consists of a six-layer all-stereo inner drift chamber and a 47-layer main drift chamber. 
These chambers are coaxial with a superconducting solenoid that provides a uniform 1.0- 
Tesla magnetic field throughout the volume occupied by all active detector components used 
for this analysis. Charged particles are required to satisfy criteria ensuring successful fits 
and vertices consistent with the e~^e~ collision point. The resulting momentum resolution 
is ~ 0.6% at 1 GeV/c for tracks that traverse all layers of the drift chamber. Oppositely- 
charged and vertex-constrained pairs of tracks are identified as Kg vr+vr" candidates if 
their invariant mass is within 4.5 standard deviations (a) of the known mass (~ 12 MeV/c^). 

The main drift chamber also provides dE/dx measurements for charged-hadron identifi- 
cation, complemented by a Ring-Imaging Cherenkov (RICH) detector covering 80% of Air. 
The rate of pions faking kaons is (1.10 ± 0.37)%, with a pion identification efficiency for 
tracks in the RICH of (94.5 ± 0.4)%. The rate of kaons faking pions is (2.47 ± 0.38)%, with 
a kaon identification efficiency for tracks in the RICH of (88.4 ± 0.6)%. 

An electromagnetic calorimeter consisting of 7784 CsI(Tl) crystals provides electron iden- 
tification and neutral detection over 93% of Att, with photon-energy resolution of 2.2% at 
1 GeV and 5% at 100 MeV. We select tt'' and 1] candidates from pairs of photons with in- 
variant masses within 3a of the known values |5j (a ~ 6 MeV/c^ for 7r° and a ~ 12 MeV/c^ 
for rj). 

TABLE I: Center-of-mass energies and integrated luminosity totals for all data samples used in 
this paper. 



Scm (MeV) 


JC dt (pb-i) 


3970 


3.85 


3990 


3.36 


4010 


5.63 


4015 


1.47 


4030 


3.01 


4060 


3.29 


4120 


2.76 


4140 


4.87 


4160 


10.16 


4170 


178.89 


4180 


5.67 


4200 


2.81 


4260 


13.11 
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TABLE II: The decay modes and branching fractions used in determining the Dg cross sections. 



Modes 


i3 (%) 


K+K-TT+, \Mkk - M^l < (10 MeV/c2) [15j 

rip^,r] 77, /9+ — > tt+tt'^ [5] 
rj'n'^ ,r]' vr+vr^?], — > 77 [SJ US] 
ry'/a"*", 7/' — > TT+vr^ry, r/ — > 77, — > 7r"'"7r'' [5] 
(pp+,(b^ K+K' , p+ ^ 7r+7r° 
KsK+,Ks vr+vr- [51 ES] 


1.99 ±0.11 

9 9 _i_ n 

z.z zn u.u 
0.62 ± 0.08 

4.3 ± 1.2 
0.66 ±0.07 

1.8 ±0.5 

3.4 ± 1.2 
1.03 ±0.06 



III. EVENT-SELECTION PROCEDURES 

The procedures and specific criteria for the selection of D^, and mesons closely 
follow previous CLEO-c analyses and are described in Refs. [12| and [15]. Candidates are 
identified based on their invariant masses and total energies, with selection criteria optimized 
on a mode-by-mode basis. We use only the cleanest final states for [K~7r^) and 
{K~7T~^7i^) selection, since these provide sufficient statistics for precise cross-section deter- 
minations. For we optimize for efficiency by selecting the eight decay modes listed in Ta- 
ble lll| Accepted intermediate-particle decay modes (mass cuts) are (p — > K^K~ (±10 MeV), 
K K-n+ (±75 MeV), r]' r/TT+Tr" (±10 MeV), and p+ ^ tt+ttO (±150 MeV). 

To determine the production yields and cross sections for the final states accessible at a 
particular center-of-mass energy, we classify events based on the energy and momentum of 
/^(s) candidate in the form of energy difference (AE = Ed^^^ — i^beam) and beam-constrained 

mass (Mbc = y^-Ebcam ~ l-^-D(^) P)- Fig. [ij shows the expected behavior in a two-dimensional 
plot of AE vs. Mbc for a Monte Carlo simulation of CLEO-c data at 4160 MeV with about 




1.80 1.90 1.95 2.00 2.10 



Mbc(GeV/c2) 

FIG. 1: AE vs. Mbc in a Monte Carlo simulation of CLEO-c data at a center-of-mass energy of 
4160 MeV, showing clear separation among the expected two-charm-meson final states. 
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FIG. 2: Momentum spectra in Monte Carlo (top) and data (bottom) at 4160 MeV for D'^ K~-k^ 
candidates with an invariant mass within 15 MeV of the nominal value. As described in the text, 
the concentrations of entries correspond to expected final states with two charm mesons. 

ten times the statistics of our data sample at that energy. There is clear separation of 
events among the expected final states consisting of two charm mesons. This separation was 
exploited during the scan run for a fast-feedback "cut-and-count" determination of event 
yields. It is also evident in plots of the momenta of charm-meson candidates selected by 
cutting on candidate invariant mass that the composition of final states can be analyzed by 
fitting the momentum spectra of , and candidates. Fig. [2] illustrates this with the 
momentum spectra for K^-k^ candidates within 15 MeV of the nominal mass both 

in the Monte Carlo sample of Fig. [l| and in 10.16 pb"^ of CLEO-c data at 4160 MeV. While 
no background corrections have been applied to these distributions, the structure of distinct 
Doppler-smeared peaks corresponding to different final states is evident. The Monte Carlo 
and data show good qualitative agreement, with concentrations of events corresponding to 
prominent final states near 0.95 GeV/c {DD), 0.73 GeV/c {D*D) and 0.5 GeV/c {D*D*). 

The cross sections for all contributing final states can be determined by correcting the 
raw measured momentum spectra like Fig. [2] for combinatoric and other backgrounds and 
then fitting to Monte Carlo predictions of the spectra. To achieve good fits, all significant 
production mechanisms must be included and the spectra predicted by Monte Carlo must 
refiect correct D*-decay angular distributions and the effects of initial state radiation (ISR). 

IV. EVIDENCE FOR MULTI-BODY PRODUCTION 

While the qualitative features of the measured charm-meson momentum spectra accorded 
with expectations (Fig. [2]), initial attempts to fit the spectra did not produce acceptable 
results. It was quickly concluded that the two-body processes listed above are insufficient to 
account for all observed charm-meson production. Final states like DD^*^7i{Tr . . .), in which 
the charm- meson pair is accompanied by one or more additional pions, emerged as the likely 
explanation. While not unexpected, these "multi-body" events have not previously been 
observed in the charm-threshold region, and there are no predictions of the cross sections 
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FIG. 3: The mass spectrum of X in (a) e+e" ^ L>°7r=^X at 4170 MeV, (b) e+e" ^ D*^tt^X 
at 4170 MeV, and (c) e+e" ^ L>*°7r=^X at 4260 MeV. Peaks at the D* mass in (a) and the D 
mass in (b) are evidence for the decay D*Dtt. The D peak in (c) confirms D*Dtt and the D* peak 
demonstrates that D*D*tt is produced at 4260 MeV. 

for and production through multi-body final states. 

To assess which multi-body final states {DDtt, D*Dtt, etc.) are measurably populated 
in our data, we examine observables other than the charm-meson momenta, because ISR 
causes smearing of the peaks in the momentum spectra that can obscure the two-body 
kinematics. We applied D^*'^ momentum selection criteria to exclude two-body contributions 
and examined the distributions of missing mass against a D^*'^ and an accompanying charged 
or neutral pion, using charge correlations to suppress incorrect combinations. Fig. [3] shows 
clear evidence for D*Dtt events at 4170 MeV, as well as indications of D*D*tt in the sample 
of 13 pb~^ collected at 4260 MeV (Fig. [sj^). These events cannot be attributed to two-body 
production with ISR, because radiative photons would destroy any peaking in the missing- 
mass spectrum. The absence of a peak at the D mass in Fig. [3^ indicates that there is 
no evidence for DDtc production. Analysis of events with Dg reveals no evidence for multi- 
body production, consistent with expectations, since the D^D~-n^ final state violates isospin 
conservation. 



-a) 




- . X n 

- b) L 

4. 




\ — ' — 1 — ' — \ — ' — \ — 1- 
_ c) 


J 


J\ - 



7 



V. MOMENTUM-SPECTRUM FITS AND CROSS-SECTION RESULTS 



Momentum spectra for D^, and candidates were found by requiring the invariant 
mass of decay products to be within ±15 MeV of the nominal value. Backgrounds 
are estimated with a sideband technique. Sideband regions are taken on both sides of the 
expected signal, and are significantly larger than the signal region to minimize statistical 
uncertainty in the background subtraction. Sideband widths are set mode by mode based 
on expectations for specific background processes. 

Having identified the components of multi-body charm production, we determine yields 
for these channels and the two-body modes by fitting the sideband-subtracted D^, and 
momentum spectra. Signal momentum distributions for specific channels are based on 
full GEANT [16] simulations using EvtGen p.7j for the production and decay of charm mesons. 
The EvtGen simulation incorporates all angular correlations by using individual amplitudes 
for each node in the decay chain. ISR is included in the simulation, which requires input of 
energy-dependent cross sections for each final state. We used simple parameterizations of 
these cross sections constructed by linearly interpolating between the preliminary measure- 
ments from our analysis. (In doing this we made the assumption that the energy dependence 
of the Born-level cross sections is adequately represented by the uncorrected cross sections.) 
For the multi-body D*Dn and D*D*7i final states we used a spin-averaged phase-space 
model within EvtGen. 

The momentum- dependent yields and fits to the relatively large sample of data at 
4170 MeV are shown in Fig.Qfor (a) D° K~7r+, (b) D+ K-n+n+, and (c) D+ 07r+ 
candidates. The lack of D'^ entries below ~ 200 MeV confirms the absence of multi-body 
Ds production. Because of the relative simplicity of Dg production, demonstrated by the 
(pTT^ fits, and the limited statistics of the sample, we determine the final cross sec- 
tions for D'^D~ , D*~^D~, and D*'^D*~ by using a sideband-subtraction technique to count 
signal events in a region of the Mbc - AE plane. The cross sections are then determined 
from a weighted sum of the yields for the eight Dg decay modes given in Table [IT], with 
weights minimizing the combined statistical and systematic uncertainties calculated from 
previously measured branching fractions and efficiencies determined by Monte Carlo. The 
cut-and-count analysis gives results that are consistent with momentum fits. There is good 
agreement among the separately-calculated cross sections for the different Dg decay modes. 

Each of the thirteen data subsamples has been analyzed with the techniques developed 
and refined on data at 4170 MeV. A complete set of fit results is provided in Ref. [T8] . 
Fig. [5] shows the D^, and Dg fits for data sample at 4260 MeV, which are of particular 
interest because the charm-production cross sections might provide insight to the nature of 
the y(4260) state. The fits at 4260 MeV behave similarly to those at lower energy, although 
a larger proportion of multi-body decays is apparent. 

Cross sections for the two-body and multi-body final states are shown in Fig. [6] The 
uncertainties on the data points are statistical and systematic combined in quadrature. 
Ref. |18] provides detailed descriptions of the systematic uncertainties of the cross-section 
determinations. Briefiy, there are three sources of systematic uncertainty: determination of 
the efficiency of charm-meson selection, extraction of yields, and overall normalization. The 
total systematic uncertainty (Table III) is not dominated by any one of these. 

Track selection and particle identification closely follow previous CLEO-c analyses [12], 
[15] . The efficiency for reconstructing charged tracks has been estimated by a missing- mass 
technique applied to events collected at the i'i^S) and '?/'(3770) resonances. There is good 
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FIG. 4: Sideband-subtracted momentum spectra for (a) K^tt^, (b) K~tt~^it'^, and 

(c) Df — > (f)7r~^ at 4170 MeV. Data are shown as points with errors and the total fit result is shown 
as the solid black line. The colored histograms represent specific -D(g)-production mechanisms, 
with shapes obtained from Monte Carlo simulations and normalizations determined by the fits. 
For example, the primary in D*^D^, which peaks at 0.7 GeV/c, is shown in bright red. The 
secondary mesons from the primary D*^ decaying via the emission of a vr'^ form the broad peak 
at 0.6 GeV/c shown in light blue. The second broad peak, at 0.6 GeV/c, consists of mesons 
from the charged pion decay of the D*+ in D*^D^ . All sources of multi-body events are combined 
and result in the broad spectrum between and 0.5 GeV/c shown in dark red. 



agreement between data and Monte Carlo, with an estimated relative uncertainty of ±0.7% 
per track. Pion and kaon identification has been studied with D° and decays in -^(3770) 
data, with estimated systematic uncertainties in the respective efficiencies of ±0.3% and 
±1.3%. The uncertainties on reconstruction efficiencies for the neutral particles 7r° and rj 
for Ds decays have been estimated at ±2% and ±4%, respectively. 
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FIG. 5: Sideband-subtracted momentum spectra for (a) K^tt^, (b) K^tt^tt^, and 

(c) Df — > (pTT^ at 4260 MeV. Data are shown as points with errors and the total fit result is 
shown as the solid black line. The colored histograms represent fit components, mostly single D(^gy 
production modes. For example, the primary in D*^D^, which peaks at 0.8 GeV/c, is shown 
in bright red. The secondary mesons from the primary D*^ decaying via the emission of a vr'^ 
form the broad peak at 0.7 GeV/c shown in light blue. The second broad peak, at 0.7 GeV/c, 
consists of mesons from the charged pion decay of the D*^ in D*~^D~ . The multi-body events 
are combined and result in the broad spectra between and 0.6 GeV/c for D*D7r (dark red) and 
between and 0.4 GeV/c for D*D*tt (black). 



The extraction of event yields by fitting the charm-meson momentum spectra (non-D^ 
modes) incurs systematic uncertainty primarily through the signal functions generated by 
Monte Carlo, which depend on details of ISR and, in the case of D*D*, the helicity ampli- 
tudes [18] and resulting D-meson angular distributions. As for the exclusive measurements, 
these details were studied with the large data sample at 4170 MeV, for which statistical 
uncertainties are small, and the resulting estimated relative systematic uncertainties are 
applied to all energy points. For the ISR calculation, the exclusive cross sections input to 
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FIG. 6: Exclusive cross sections for two-body and multi-body charm-meson final states, and total 
observed charm cross section with combined statistical and systematic uncertainties. 



EvtGen were varied from their nominal shapes. While a qualitative constraint of consistency 
with our measured cross sections was imposed, some extreme variations are included in the 
final systematic uncertainty. Both the direct effect on the fitted yield of varying a specific 
mode and the indirect effect of varying other modes were computed, although the former 
dominates in quadrature. 

The yields for Dg final states are determined by direct counts after cutting on Mbc and 
AE. Systematic uncertainty arises in these measurements if the Monte Carlo simulation 
does not provide an accurate determination of the associated efficiency. This is probed by 
adjusting the selection criteria and recomputing the cross sections, again using the high- 
statistics sample at 4170 MeV. The systematic uncertainties assigned based on these studies 
are ±3%, ±2.5% and ±5% for D+D~, D*+D~, and D*+D*', respectively. 

In converting the measured yields to cross sections we must correct for the branching 
fractions of the charm-meson decay modes. For each of the non-strange charm mesons. 
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TABLE III: Total systematic errors on the exclusive cross sections. 



A/Tnrlp 


IVcidLlVc J_jliui \ -Lvy } 


1 Iti'f" tiT'TTi m ti/i IrwT ^/l r~»TTi tiirt" n in H Tt"G 
X^c Icl llilllcvi uy iVlUllicllL Ulli 1 1 Lo 




JJU 


4.0 


n n* 

UU 




D* h* 


4 7 


D*Dtt 


12.0 


D*D*TT 


25.0 


Determined by Counting 






5.6 




5.3 


Dl+D*- 


6.8 



only one mode is used and CLEO-c measurements [12] provide the branching fractions and 
uncertainties: ±3.1% for K^tt^ and ±3.9% for K^ti^ti^. For Dg modes we 

use CLEO-c measurements of the branching fractions for the eight decay modes included 
in the weighted sum |l5j. The world-average value is used for the D*^ D^tt~^ branching 
fraction, with a systematic uncertainty of ±0.7% [5]. Finally, the cross-section normalization 
also depends on the the absolute determination of the integrated luminosity for each data 
sample, with a systematic uncertainty of ±1.0% [13]. 

A mode-by-mode summary of the systematic uncertainties in the exclusive cross-section 
measurements is provided in Table III, The systematic errors are 100% correlated across 
energy. 

The cross-section measurements are presented in Tables IV, V, VI (modes with only two 
charm mesons), and VII (multi-body modes). 

As a cross-check, for the two largest data samples (4170 MeV and 4260 MeV), the multi- 
body cross sections are also determined by fitting the distributions of missing mass against 
detected D^n, D^n and D*7i combinations. While these measurements are less precise, they 
show good agreement with the results of the momentum-spectrum fits. 



VI. INCLUSIVE CROSS-SECTION MEASUREMENTS 

If all final states have been included, the sum of the exclusive cross sections should equal 
the total charm cross section. We test this supposition with two inclusive measurements 
that can also be compared with past results. 

The first cross-check is a measurement of the total charm-meson cross section: 

^(e+e ^DDX) = ^, (1) 

where the contributing cross sections are defined by cxd = N^/eBC, where e and B are the 
efficiency and branching fraction for the decay mode used {D^ K'n'^, — > K~TT~^'n'~^, 
and K~ K~^n^), L is the integrated luminosity, and is the yield obtained by fit- 

ting the mass spectrum. In the case of and D"'", the invariant-mass distribution is fitted 
to a Gaussian signal and polynomial background. For D^, the event-type requirements are 
maintained because of the relatively large background for the high-yield K^K^-n^ decay 
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TABLE IV: Measured cross sections for final states consisting of two neutral non-strange charm 
mesons. The first error on each cross section is statistical and the second is systematic. 





(J [J-^ J-' J 


/7-)*0 p)0\ 
(J I / y 1 J J 
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I IVlc V J 








oy / u 


oo in zy in ^ 


ooen -1- 1 "iA -1- 78 
zzou in lo^ m / o 






TOO _1_ A-l _l_ a 

loo in ^1 HZ u 


z 1 4:U HI lo i m yo 








oozu HI lo m iio 








oo^u HI zoo HI i-Oi 


ZIO HI / U HI y 


4030 


iJiJt: —1— 1 \J —1— 


Qonn +183+109 


onnn + i 05 + Q4 

^KjyjKJ _1_ i^^fj _1_ t/T: 


4060 


410 ± 72 ± 18 


2230 ± 147 ± 76 


2290 + 132 + 108 


4120 


303 ± 70 ± 14 


1400 ± 135 ± 48 


2550 + 154 + 120 


4140 


177 ± 40 ± 8 


1350 + 100 + 46 


2443+ 116+ 115 


4160 


167 ± 28 ± 8 


1252 + 69 + 43 


2566 + 84+ 121 


4170 


177±7±8 


1272 + 19 + 43 


2363 + 19+ 111 


4180 


179 ± 39 ± 8 


1211 + 92 + 41 


2173+ 104+ 102 


4200 


180 ± 55 ± 8 


1030 + 123 + 35 


1830 + 139 + 86 


4260 


86± 18±4 


1080 + 59 + 37 


269 + 42+ 13 



TABLE V: Measured cross sections for final states consisting of two charged non-strange charm 
mesons. The first error on each cross section is statistical and the second is systematic. 





a{D+D-) 


a{D*+D-) 


(j{D*+D*-) 


(MeV) 


(Pb) 


(pb) 


(pb) 


3970 


137 + 26 + 6 


2230 + 131 + 76 




3990 


90 + 22 + 4 


2750 + 156 + 94 




4010 


135 + 22 + 6 


3300 + 132 + 112 




4015 


38 + 20 + 2 


3703 + 274 + 126 




4030 


196 + 35 + 9 


3300 + 181 + 112 


1400 + 170 + 66 


4060 


480 + 55 + 22 


2170 + 143 + 74 


2390 + 222 + 112 


4120 


310 + 50+14 


1560 + 136 + 53 


2280 + 232 + 107 


4140 


200 + 29 + 9 


1376 + 98 + 47 


2556 + 196 + 120 


4160 


200 + 21 + 9 


1376 + 69 + 47 


2479 + 135 + 117 


4170 


182 + 6 + 8 


1285+ 18 + 44 


2357+ 19 + 111 


4180 


197 + 27 + 9 


1296 + 87 + 44 


2145 + 172 + 101 


4200 


181 + 36 + 8 


1070 + 116 + 36 


1564 + 215 + 74 


4260 


94+13 + 4 


1022 + 54 + 35 


237 + 54+ 11 



mode. For our energy points below 4120 MeV, where Dg production occurs only through 
DfDj, the yield is extracted by fitting Mbc to a Gaussian signal and ARGUS background 
function [19]. For 4120 MeV and above, event types involving D*~^ contribute. For all can- 
didate events that pass the selection requirements for any of D'^Dj, D*^Dj , and D*^D*~ 
(the last only for 4260 MeV), a fit to the invariant mass is used to determine the yield. 



13 



TABLE VI: Measured cross sections for final states consisting of two strange charm mesons. The 
first error on each cross section is statistical and the second is systematic. 



-'-'era 


rr( n+ n^) 




"\^s ) 


I IVlc V J 


IP'-'; 






oy / u 










1 QQ -1- Q1 -1-7 
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4030 


-Lit: —1— KJ\J —1— VJ 






4060 


51 ±28 ±3 






4120 


26 ± 26 ± 1 


478 ± 64 ± 25 




4140 


25 ± 20 ± 1 


684 ± 59 ± 36 




4160 


< 15 (90%C.L.) 


905 ± 11 ±48 




4170 


34 ± 3 ± 2 


916 ± 11 ±49 




4180 


7±16±1 


889 ± 59 ± 47 




4200 


15 ± 22 ± 1 


812 ±82 ±43 




4260 


47 ± 22 ± 3 


34 ± 9 ± 2 


440 ± 27 ± 30 



TABLE VIL Measured cross sections for multi-body final states, consisting of two charm mesons 
and an extra pion, for all data points above the production threshold. The first error on each cross 
section is statistical and the second is systematic. 



^cm (MeV) 


a{D*D7:) (pb) 


a{D*D*7r) (pb) 


4060 


144 ± 94 ± 17 




4120 


45 ± 83 ± 5 




4140 


412 ± 87 ± 49 




4160 


389 ± 60 ± 47 




4170 


440 ± 20 ± 53 




4180 


575 ± 92 ± 69 




4200 


735 ± 129 ± 88 




4260 


638 ± 93 ± 77 


322 ±67 ±80 



The second cross-check is a determination of the total cross section made by counting 
multihadronic events. The contribution of uds continuum production is estimated with 
measurements made at E^ui = 3671 MeV, below cc threshold, and extrapolated as 1/s. 
Procedures for this measurement are identical to those used to determine the cross section 
for e+e- ^/'(3770) ^ hadrons in CLEO-c data at E^^a = 3770 MeV [20]. 

Fig. [6] (bottom frame) shows the inclusive measurements (statistical and systematic un- 
certainties combined in quadrature) and the sum of the cross sections for the measured 
exclusive final states without radiative corrections. The excellent agreement demonstrates 
that, to current precision, the measured exclusive two- and three-body final states saturate 
charm production in this region. Furthermore, charm is demonstrated to account for all 
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production of multihadronic events above the extrapolated uds cross section. 

For the inclusive-charm cross-section measurements, the systematic uncertainties associ- 
ated with the per-particle efficiencies for tracking and particle identification are identical to 
those of the exclusive measurements. The uncertainties in normalization (luminosity and 
branching fractions) are also identical. Systematic uncertainty in the yield extraction is 
dominated by the choice of fitting function. This is evaluated mode by mode and prop- 
agated into overall systematic uncertainties accounting for all correlations, with combined 
systematic uncertainties of ±4.3%, ±5.1%, and ±8.6% (±10.6%) for D^, D+, and D+ below 
(above) 4120 MeV. For the hadron-counting inclusive cross sections, the systematic uncer- 
tainties are identical to those of Ref. [20j. The systematic errors for the hadron-counting 
method are slightly energy-dependent, varying between 5.2% and 6.1% due to the different 
amounts of J/ip, ip{2S), and ?/'(3770) present at each energy. 



Table VIII gives the inclusive cross sections and the sum of the exclusive cross sections 



with both statistical and systematic uncertainties. 

TABLE VIII: Comparison of the total charm cross section determined by summing the exclu- 



sive measurements (Tables IV VI and VII) with those found by the two inclusive techniques: 
charm-meson counting and multihadronic-event counting. The first error on each measurement 
is statistical and the second systematic. The cross-section measurements are not radiatively cor- 
rected. The last column gives the value of R from the hadron-counting measurement, with radiative 
corrections as described in the text and correction for non-charm continuum production based on 
Ruds = 2.285 ± 0.03, as determined by a - fit to previous R measurements between 3.2 and 3.72 
GeV L23J. 



Energy 


Exclusive 


Inclusive 


Hadron 




R 




(MeV) 


D-meson (nb) 


L)-meson (nb) 


Counting (nb) 


(ISR-corrected) 


3970 


4.83 ±0.19 ±0.15 


4.91 ±0.18 ±0.16 


4.91 ±0.13 ±0.30 


3 


36 ± 0.04 ±0 


05 


3990 


5.85 ±0.23 ±0.19 


5.93 ±0.21 ±0.19 


5.87 ±0.14 ±0.34 


3 


55 ± 0.05 ±0 


06 


4010 


7.10 ±0.14 ±0.23 


7.05 ±0.17 ±0.23 


7.21 ±0.12 ±0.40 


3 


88 ± 0.04 ± 


08 


4015 


7.94 ±0.41 ±0.26 


7.62 ±0.34 ±0.25 


7.88 ±0.18 ±0.43 


3 


95 ± 0.08 ±0 


08 


4030 


10.60 ±0.34 ±0.27 


10.87 ±0.28 ±0.37 


11.30 ±0.15 ±0.59 


4 


74 ± 0.07 ± 


12 


4060 


10.16 ±0.36 ±0.27 


9.98 ±0.26 ±0.34 


9.98 ±0.14 ±0.53 


4 


34 ± 0.05 ± 


10 


4120 


8.95 ±0.37 ±0.25 


9.13 ±0.28 ±0.31 


9.43 ±0.15 ±0.49 


4 


21 ± 0.06 ±0 


10 


4140 


9.22 ±0.29 ±0.26 


9.11 ±0.22 ±0.30 


9.58 ±0.24 ±0.50 


4 


18 ± 0.04 ±0 


10 


4160 


9.33 ±0.20 ±0.26 


9.10 ±0.15 ±0.30 


9.62 ±0.17 ±0.50 


4 


18 ± 0.03 ±0 


10 


4170 


9.03 ±0.04 ±0.25 


9.09 ±0.07 ±0.30 


9.45 ±0.09 ±0.49 


4 


20 ± 0.01 ±0 


10 


4180 


8.67 ±0.27 ±0.24 


8.70 ±0.20 ±0.29 


9.07 ±0.12 ±0.47 


4 


17 ± 0.04 ±0 


10 


4200 


7.42 ±0.35 ±0.20 


7.45 ±0.26 ±0.25 


8.37 ±0.14 ±0.43 


3 


77 ± 0.05 ±0 


08 


4260 


4.27 ±0.16 ±0.14 


4.20 ±0.10 ±0.14 


4.34 ±0.16 ±0.23 


3 


06 ± 0.02 ±0 


04 



For comparison with other experiments and theory it is necessary to obtain Born-level 
cross sections from the observed cross sections by correcting for ISR. We do this by cal- 
culating correction factors following the method of Kuraev and Fadin [21] , which gives the 
observed cross section at any y/s: 

1 

<7obs(s) = / dk ■ f{k,s)aB{ses), (2) 
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FIG. 7: R (including radiative corrections) from this analysis and from previous measurements 

where the Bom cross section is a function of the effective center-of-mass energy squared 
{k = (s — Ses)/s), and f{k,s) is the ISR kernel. The radiative-correction factor is also 
calculated following the alternative implementation of Bonneau and Martin [22j . We take the 
difference between the two methods as an estimate of the underlying theoretical uncertainty 
in the calculation of the radiative-correction factor. We also consider systematic uncertainty 
due to our approximation of crB{ses), required for Eq. 2, by taking the difference between 
a simple linear interpolation and a fit to a sum of Breit-Wigners to both the BES [6J and 
Crystal Ball (CB) [23j R measurements. Fig. [T] shows that there is excellent agreement 
between our inclusive-charm measurement and the previous R measurements. 

VII. SUMMARY AND CONCLUSIONS 

In summary, we have presented detailed information about charm production above cc 
threshold. Realizing the main objective of the CLEO-c scan run, we find the center-of- 
mass energy that maximizes the yield of Dg to be 4170 MeV, where the cross section of 
~ 0.9 nb is dominantly D*^D~ . This information has guided the planning of subsequent 
CLEO-c running, with initial results already presented on leptonic [21] and hadronic [I5] Ds 
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FIG. 8: Comparisons between measured cross sections and the updated predictions of the potential 
model of Eichten et al. P, f25| (solid lines). 



decays. The total charm cross section between 3.97 GeV and 4.26 GeV has been measured 
both inclusively and for specific two-body and multi-body final states. Internal consistency 
is excellent and radiatively-corrected inclusive cross sections are consistent with previous 
experimental results. Fig. [g] shows that the observed exclusive cross sections for DD, D*D, 
D*D*, DfD-, D*+D-, D*+D*-, D*D7i, and D*D*7i exhibit structure that reflects the 
intricate behavior expected in the charm-threshold region. Fig. |8] provides a comparison 
between our measured cross sections and the updated calculation of Eichten et al. [9l |25] . 
There is reasonable qualitative agreement for most of the two-charm-meson final states. The 
most notable exception is the cross section for D*D* in the region between 4050 and 4200 
MeV, where the measurement exceeds the prediction by as much as 2 nb. This corresponds 
to nearly a factor-of-two disagreement in the ratio of D*D* to D*D production, accounting 
for about two thirds of the difference in the total charm cross section. This is a much larger 
effect than the absence of a multi-body component from the theoretical prediction. 

It has been suggested by Dubynskiy and Voloshin [26] that the existence of a peak in the 
D*D and D^D~ channels at the D*D* threshold, along with the observation that there is a 
minimum in DD, in agreement with recent results from BaBar [27], can be interpreted as a 
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possible new narrow resonance, but available data are insufficient for a definitive assessment. 

The D*D* cross section exhibits a plateau just above its threshold. This contrasts with 
D*D, which we observe to peak at threshold, in agreement with recent results from Belle 



Studies of open-charm production at 4260 MeV have the potential to discriminate among 
possible explanations of the nature of the y(4260). For example, hybrid charmonium mod- 
els predict a large coupling to the wide 1)1(2430)"/)° and a small one to D^D~ |29]. A 
tetraquark interpretation suggests a large decay to DD or D^D~ [291 EHl EI] ■ Complicated 
threshold effects could lead to enhancement of the D*ti final state through off-shell pro- 
duction of Di [32]. Tables IV, VI and VII show no evidence for enhancement of the cross 



section for any open-charm final states at 4260 MeV. CLEO-c has previously confirmed the 
F(4260) through its decay to tt+tt" J/-?/;, measuring a{7i~^7i~ J/ip) = 58t\l ± 4 pb [2J. Under 
the assumption that all open-charm production is accounted for by y(4260) decays, it is 
possible to set conservative upper limits on the ratio of the cross section for production of 
y(4260) and decay to our measured open-charm states to that for production and decay to 
n~^n~J/il). Table IX provides a compilation of these limits. The lack of obvious enhance- 



TABLE IX: Upper limits (90% confidence level) on the ratio of the cross section for production 
of y(4260) and decay to our measured open-charm states at 4260 MeV to that for production of 
y(4260) and decay to n+ir-J/i;. 



Final State (X) 


cr(y(4260)^X) 


o-(y (4260)^7r+7r- J/^/)) 


DD 


< 4.0 


D*D 


< 45 


D*D* 


< 11 


D*Dtt 


< 15 


D*D*TT 


< 8.2 


DtD- 


< 1.3 




< 0.8 


D*+D*- 


< 9.5 



ment in any open-charm channel relative to other energies, which is in stark contrast to the 
clear enhancement in vr+vr" J/?/', tends to disfavor the hybrid charmonium and tetraquark 
proposals. More definitive statements will require additional data from future experiments. 
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